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Abstract. - We study the sympathetic cooling of a trapped Fermi gas interacting with an 
ideal Bose gas below the critical temperature of the Bose-Einstein condensation. We derive 
the quantum master equation, which describes the dynamics of the fermionic component, and 
postulating the thermal distribution for both gases we calculate analytically the rate at which 
fermions are cooled by the bosonic atoms. The particle losses constitute an important source 
of heating of the degenerate Fermi gas. We evaluate the rate of loss-induced heating and derive 
analytical results for the final temperature of fermions, which is limited in the presence of 
particle losses. 



Evaporative cooling has proven to be an essential tool to obtain degenerate Fermi gases 
[1-6]. In Fermi systems, the s-wave collisions of indistinguishable particles are forbidden 
due to the antisymmetry requirement, and therefore the only way to cool fermions using 
collisions, is to cool them sympathetically by bringing them in contact with atoms in different 
hyperfine states, or another species. In the regime of quantum degeneracy, however, also the 
collisional processes in Fermi-Fermi mixtures are strongly suppressed due to Pauli blocking. As 
a consequence, the efficiency of sympathetic cooling is reduced, and the particle losses, which 
provide an important source of heating in the degenerate regime [7] , prevent from reaching the 
lower temperatures T, necessary, for instance, for achievement of the superfluid BCS state. 
This problem may be avoided by adiabatic crossing of a Feshbach resonance, since very low 
temperatures may be achieved [8] when moving from a molecular BEC regime [9] to the BCS 
regime. In this crossover, which has recently attracted a large attention, strong evidences 
for superfluidity have been observed [10-14]. Another obstacle appears when a Fermionic 
gas is sympathetically cooled using a Bose-Einstein condensate, since as a consequence of the 
superfluid properties of the condensate, the sympathetic cooling is limited to velocities larger 
than the sound velocity [15]. 

Sympathetic cooling has been recently considered by several theoretical groups [15-22]. It 
should be stressed, however, that the limitations due to particle losses have been addressed 
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in detail only recently by Carr et al. [23], who have studied the cooling of fermions by an 
ideal homogeneous Bose gas at T = 0. In this paper we extend the result of Ref. [23], and 
investigate the influence of the particle losses on the sympathetic cooling of harmonically 
trapped gases at finite T. We consider here a simplified model which assumes that fermions 
are in contact with an ideal Bose gas, which is kept at constant T . This can be achieved if 
the bosons are continuously cooled by means of evaporative cooling or eventually when the 
bosonic component is much larger than the fermionic one. 

The gas of Nf fermions is confined in a harmonic trap of frequency ojf. In the following 
we assume that the trapping potential is spherically symmetric, however, our results can 
be generalized to anisotropic traps. The fermions interact with a gas of TV^ bosonic atoms, 
confined in a harmonic trap of frequency w;,. We treat the gas of bosons as a reservoir, which 
assumes that its temperature and statistical properties do not change during the cooling 
process. The Hamiltonian of the system reads 

= ^ enflnfln + ^ em^L^m + ^ C/n,n',m,m' ajjAn' 6j„6m' , (1) 
n m n,n',m,m' 

where n and m label the states of the three-dimensional harmonic oscillator, £„ and 
denote single particle energies of fermions and bosons, respectively, and On (6m) is the creation 
operator of fermions (bosons) fulfilling usual anticommutation (commutation) relations. In 
the low-temperature regime considered here, the scattering between atoms of both species is 
purely of s-wave type and can be modeled by a delta pseudopotential with scattering length 
a. Thus the matrix elements of the boson-fermion interactions are given by 

C/„,„',m,m' j dVV';(r)i/;„'(r)0;,(r)(/)m'(r), (2) 

where '!/'n(r) and 0m (r) are the single-particle eigenfunctions of fermions and bosons, respec- 
tively, g = 2Trh?a/ fijrn and /i^ = mijirif / {mi, + 'mf) denotes the reduced mass {mi, is the mass 
of boson and mf is the mass of fermion). We apply standard techniques of the theory of 
quantum-stochastic processes [24, 25] and obtain the quantum Master equation (ME) for the 
density matrix p of the fermions [16] 

^ = [Hp + Hp_p,p]+ r^^'f i^al^a^, pal^a^' - al^a^i al^a^, p - pal^a^' al^a^,) , 

n,n' ,m,m' 

(3) 

where Hp — eninOn and the rate coefficients are of the form: 

^Z'^ = ■^'5(£n+em,en'+emOXl^"^"'^'l^kC/m,m'.k.q«q(l+4)'5(em+4:em'+£q). (4) 

In Eq. S{e,e') denotes the Kronecker delta function, which accounts for the conservation 
of energy during collisions between fermions and bosons, and denote the occupation num- 
bers of the bosonic reservoir. The operator Hp^p describes an effective interaction between 
fermions, which is mediated by bosonic atoms: 

Hf-f^ A™„? a]^a„/a]„am', (5) 

n,n',m,m' 

where 

Am.m' <-/ , , \V^j-r tt , 6 n 1 ~ ^(^m + £k , £m' + £q) 
^n.n' ='j(en+£m,en'+em')2^'^n,n',q,kt'm,m',k,q"q(l+»^k) ^ _ —f, • (6) 
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We assume that interactions between bosons and fermions are sufficiently small, and hence 
Hp-F does not influence signiflcantly the single-particle eigenstates in a harmonic trap. 

To calculate the dynamics of a fermionic gas from the ME we introduce some further 
approximations. First, by considering the decoherence induced by the interactions between 
bosons and fermions one can argue that the nondiagonal elements of the density matrix in the 
energy representation decay on a time-scale much shorter than the time scale of the dynamics 
due to sympathetic cooling [20]. Inclusion of Hp-p lifts the degeneracy of the states of 
equal energy but different occupation numbers of single-particle eigenstates, and hence the 
decoherence affects also the nondiagonal terms with respect to the occupation numbers. Then, 
we apply ergodic approximation assuming that for a given energy all diagonal elements of p 
are equal. This requires that the process of equilibration within a single energy shell is much 
faster than for the whole system. Thus, the density matrix p can be written in the form [20] 

^(*) = E T{N%) ^''"^'^ • 

Here, A enumerates microstates {n-n}\ of the fermionic component, and {n-a}\) represents 
the state with energy E and a given distribution of fermions. The probability that the system 
of fermions has energy E is denoted by psit), while T{Nf,E) is the number of microstates 
with energy E and number of particles Nf. Substituting p{t) given by Eq. into the ME Q, 
and neglecting the terms due to the Hf-f, we obtain the equations determining the dynamics 
of the probabilities pE{t): 

^ - 2 ^ {pF+e„-e^{t) - - PE{t) (n„(l - n^)) ^) , (8) 

n,l 

where (rin(l — n^)) ^ stands for the microcanonical average over microstates with energy E, 
and rin = aJjCn- Probabilities PE{t) determine all the thermodynamic quantities of the gas 
of fermions. In particular, its mean energy is given by {E) = EpE{t). Combining this 
expression with Eq. |(HJ), we calculate the rate F = — ^ (E) of cooling of the fermions by the 
bosonic reservoir: 

r = 2 ^ F™;ii (£„ - s^)J2PE{t) (n„(l - n^))^ . (9) 

n,m E 

Since we are mainly interested in the regime where fermions are sufficiently cold and the effects 
of heating due to the losses become important, we assume that the gas of fermions is close 
to equilibrium and evaluate PE{t) from the canonical or the grand-canonical ensemble. We 
chose the latter possibility and postulate the grand-canonical distribution for both fermionic 
component and bosonic reservoir. Substituting Eq. for the rates F^^^ into Eq. we 
arrive at 

F = E |t^m,n,k,q|^ rin(l - ?^m)r^q(l + n^) (en - em) (5(en + eq,£m + ek), (10) 

n,m,k,q 

where n„ = ze-"'- / {1 + ze-^'-), n\ = z^e' ^< / {I - z^e" , (3 = l/(fcsr), 13b = l/{kBn), 
and z, Zf, denote the fugacities of fermions and bosons respectively. First we perform sum- 
mation over degeneracies of the single-particle states. To this end we utilize the following 
approximate result 

irr |2x/ \u \xt \xt \ D{mm{e3,ei))g^ 

} , |L/m.n,k,q| (Cn , ei)(3 (£m , 62 )0 (Cq, 63)6 (Ck, 64) « .424 ' (H) 

n,m,k,q J 
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where Ub — \JhlmhiOh, aj — \JTifmfLOf^ D{e) denote the degeneracy of the energy sheU e, 
and min(e3,e4) represents the minimum of 63 and 64. This formula can be derived by tak- 
ing the semiclassical hmit of Eq. H10|l and comparing with the semiclassical quantum Boltz- 
mann equation in ergodic approximation [21,26,27]. The summation formula is valid for 
mm{e3, €4) /{hiOb) < A^min(ei, e2)/(fiw/) and min(e3, e4)/(?ki;b) < A^^min(ei, e2)/(fiw/) where 
A = aj /ab- For values of A which are not much larger or much smaller than one, this condition 
is fulfilled for collisions at sufhciently low temperatures, which take place close to the Fermi 
surface for fermions and close to the ground state for bosons. We have confirmed numerically 
the validity of formula Hll|) . and verified that it is quite accurate even for the lowest states of 
the harmonic potential. Substituting into Eq. (|10() yields 

^^2AuJb_ ne,il~ne.,)nl^il + nlj{ei-e2)D{inm{e3,ei))5{ei+e3,e2 + ei), (12) 

61,62,63,64 

where rie, = ze-f^""^ / {1 + ze-f^""^) and n^^ = Zbe-'^"''^ / {I- Zbe-'^"''^), and A = mla'^ al / {^l^af) . 
For a sufficiently large system, the summation over ei, 62, 63, 64 can be replaced by an integra- 
tion, which corresponds to taking the thermodynamic limit. This procedure is correct with 
respect to the collisions with the thermal part of the bosonic component. In the case of a 
Bose-condensed reservoir, the single-particle ground state requires a separate treatment. Let 
us first consider the rate Fth due to the interactions with the thermal atoms in the reservoir. 
Replacing sums by integrals, after some algebra we obtain 

which is valid for temperatures of fermions much smaller than the Fermi temperature Tp, 
and under assumption that the reservoir is Bose-condensed. Here C(s, a) denotes the Hurwitz 
zeta function [28]. We note that Fth does not depend on the number of fermions. This is 
is related with the fact that at low temperatures the collisions occur only close to the Fermi 
surface. On the contrary, taking into account that the number of the thermal atoms in a 
trapped Bose-condensed gas is A'^th — Ci'^) {f^ sTb / {huJb))^ , we observe that Fth is proportional 
to A'th- When the temperature of both components is equal, we expect that there will be no 
energy exchange between bosons and fermions and, as can be easily verified, this feature is 
properly taken into account by Eq. H13|l . To calculate the rate Fg describing the cooling due 
to collisions with the condensed part of the reservoir, we replace sums by integrals in Eq. (|10|l . 
putting 63 = or 64 = 0. This yields 

We note that at low temperatures the energy exchange with the condensate is proportional 
to the number of condensed atoms Nq, but does not depend on the number of fermions. 

We turn now to the description of the heating due to particle losses in the degenerate 
Fermi gas. The losses create holes in the Fermi sea, which after subsequent thermalization 
increase the temperature of the system [7]. We assume that particles are removed with 
a rate 7, which is the same for all single-particles states: {nn) = —jinn), where (rin) 
denotes the mean occupation number of the state n. For background losses, the rate 7 
does not depend on the density of fermions rif, while in the case of two- and three-body 
losses it is proportional to Uf and n'j, respectively. The removal process does not change 
the mean energy per particle: {{E)/Nf){t) = const. At temperatures T/Tp <^ 1, {E)/Nf = 
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iEF{Nf)lA+ {Tr^/2)k%T^/EF{Nf) with the Fermi energy EpiNf) w %(6iV/)i/3_ ^j^g 
creation of holes, can be equivalently described as a process that heats the system, but, in 
first approximation, does not change the number of particles. We obtain that the rate of 
heating is given by 

Th = Nf^Ep [i - ^ [T/Tpf] , (15) 

where we neglect terms of order higher than (T/Tp)'^. 

In the absence of particle losses, the Fermi component is cooled down until it reaches the 
temperature of the bosonic reservoir T^. The heating produced by the hole creation, however, 
prevents from reaching T^, and sets a constraint on the lowest temperature achievable with 
the sympathetic colling. The final temperature of the Fermi gas Tgn, at which the cooling 
is balanced by the hole heating, is determined by the condition: Fth + Fq = Fh. When the 
losses are sufficiently small, the final temperature of fermions Tgn only slightly differs from 
Tf), and in this case its value can be estimated analytically. To this end, we perform a Taylor 
expansion of the cooling rate F = Fth + Fq at T = Tb , keeping only the lowest order term 
F(r) « {dT/{dT))\T=niT - n). This yields 



12AL0bikBTb)'^T~Tb 



hujh 



(16) 



which is valid for |T — Tb\ <C Tj,. Here c is a numerical constant given by c = X]/^o C(3, 1 + 
/)C(4, 1 + ~ 1.3196. For temperatures of the bosonic reservoir Tb much smaller than the 
critical temperature Tc, one can neglect c{kBTb/ {fuvb))^ in comparison to ({3)No, and in this 
regime Eq. (|15|) reduces to 



12Ahwt (TbYT-Tb^. 



, , ' , ^A^o'. (17) 

This formula combined with the rate Fh, which at T <g; Tp can be approximated by Fh ~ 
Nf^Ep/A, leads to the following result for the final temperature of the Fermi gas 




rftn = rji + :— --f-^ . (18) 



The lowest temperature in the sympathetic cooling of a Fermi gas is achieved when the 
reservoir is totally condensed. In this particular case, the final temperature can be also 
calculated analytically. We observe that at Tb = 0, Fth = 0, while for T <^Tp the rate of loss- 
induced heating can be approximated by Fh ~ NfjEp/A. This yields the final temperature 
of fermions: 

1/3 



\ALUbNo 



Tfin = bTp I —i— I , (19) 



where b = {G^/^n/ {96({3))y/^ ~ 0.367. We note that Tfin very weakly depends on the number 
of fermions Nf, and is mainly determined by the efficiency of the cooling process (through 
parameters A'o, A, cob), and by the loss rate 7. 

To estimate the influence of the loss-induced heating on the final temperature of fermions, 
we have performed numerical calculations for the mixture of bosonic ^Li with fermionic ^Li 
[2,3]. In our calculations we assume: Nb — 1.5 x 10'', Nf — A x 10^, ujb — ujf — 27: x 1270IIz, 
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Fig. 1 Fig. 2 

Fig. 1 - Cooling rates F (solid lines) and heating rates Fh (dashed lines) versus the fermion temper- 
ature T. The rates were calculated for three different temperatures of the bosonic reservoir: Tb — 0, 
O.OSTf, 0.04rF, and for three different loss rates 7 = 0.04Hz, O.OlHz, and O.OOlHz. The values of the 
other parameters correspond to the experiment described in [2] (see text for details). 

Fig. 2 - Final temperature of the Fermi gas versus the temperature of the bosonic reservoir Tj 
calculated for the parameters of Fig. The plot shows the numerical values evaluated from the rates 
given by Eqs. (1131 . (I14II . and 1151 (solid lines), the approximate result of Eq. (1181 (dashed lines), and 
the line T — Tt, which is the final temperature in the absence of losses (dotted line). 



and a = 2.0nin, which corresponds to the experimental parameters of Ref. [2]. Fig.^shows, as 
a function of the temperature of the fermions, the coohng rates F = Fq + Fth evaluated from 
Eqs. (fl^ and (fTI|l for three different temperatures of the bosons: Tb — 0, O.OSTf, and 0.04Tf 
(solid lines). The dashed lines show the heating rate Fh given by Eq. (|15|l . Since the loss rate 
7 depends on the particular experimental conditions, we have used three different values of 7: 
0.04Hz, O.OlHz and O.OOlHz. The final temperature of fermions is determined by the crossing 
of F and Fh. We observe that for the totally condensed reservoir (T^ = 0) and for fermion 
temperatures T < O.OlTc, the cooling rates become very small due to Pauli blocking. In 
this regime the presence of even very small losses prevents from reaching lower temperatures. 
For higher temperatures of the bosonic component {Ti, = O.OSTp and = 0.04rp), the 
hole heating shifts only slightly the final fermion temperature with respect to Tb, and Tfin 
can be determined by expanding the cooling rate around Tb as in the case of Eq. ()18|) . The 
dependence of the final fermion temperature Tfin on Tb is shown in Fig. 12 where the results 
calculated from the condition Fth + Eq = Fh, using Eqs. (|13|l . (|14|l and (|15l) . are compared 
with Eq. ((TH|l . 

In our analysis we have neglected the interactions within the bosonic reservoir. This can be 
a good approximation for the conditions of experiment [2] , where the condensate is relatively 
small {Nb ^ 10"*) and the scattering length for ^Li atoms in the state \F = l,mF = —1) is 
a = 0.27nm. In general, however, the interactions modify the density profile of the bosonic 
cloud, and introduce the quantum depletion of the condensate. For repulsive forces the size of 
an interacting condensate is much larger in comparison to the ideal-gas. This can reduce the 
collisional rate between bosons and fermions and consequently decrease the cooling efficiency. 
On the other hand, the quantum depletion reduces the number of condensed atoms, transfering 
them into excited states, which increases the effective temperature of the reservoir. Therefore, 
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the most effective cooling is achieved for noninteracting bosons, which can be eventually 
realized with the help of Feshbach resonances. 

Summarizing, we have investigated the sympathetic cooling of a trapped Fermi gas in 
contact with a Bose-condensed reservoir, and shown that the loss-induced heating limits the 
cooling at low temperatures. In the regime of weak heating and in the case of totally condensed 
Bose gas, we have derived the analytical results for the final temperature of the fermions. 
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